Planar approximation for the frequencies of spin transfer oscillators 
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A large class of spin transfer oscillators use the free layer with a strong easy plane anisotropy, 
which forces its magnetization to move close to the plane. We show that in this situation the effective 
planar approximation provides a fast and accurate way of calculating the oscillator frequency. 

PACS numbers: 75.76.+j, 75.78.-n, 85.75.-d 



I. INTRODUCTION 

Spin transfer devices can have regimes in which 
their magnetic moments perform perpetual precessional 
motions /i'^ In this case they are also called spin torque 
oscillators (STO). Magnetic oscillations induced by di- 
rect current are intensively studied experimentally^"— 
and theoreticallyJ^""i2^ In the STO regime the energy is 
constantly supplied to the device from the current source 
through the spin transfer mechanism. At the same time 
it is lost through the usual dissipation mechanisms, ac- 
counted for by the Gilbert damping constant. In a state 
of steady precession the energy gain and loss are balanced 
on average. In the limit of small damping one observes 
the following general picture of the STO operation, The 
magnetic moment moves close to the trajectory which it 
would follow in the absence of damping and spin trans- 
fer. The actual trajectory is a perturbation of the zero- 
damping trajectory, chosen so as to balance the small 
dissipation with the equally small energy gain. The main 
difficulty in describing the precession states is the lack of 
knowledge about the unperturbed trajectory which is a 
solution of the complicated non-linear Landau-Lifshitz- 
Gilbert (LLG) equation. Unless one considers a small 
radius precession near an equilibrium point, the analytic 
form of such a trajectory is usually unknown and one 
is forced to use numeric methods. In this paper we will 
consider a special class of spin transfer devices with domi- 
nating easy plane anisotropy. This anisotropy often arises 
from the thin disk shape of the magnetic layers found in 
the majority of experimental structures. Due to the dom- 
inating easy plane anisotropy the LLG equation can be 
approximated by an effective planar equation^^"— that 
is less complex and easier to treat analytically. Here we 
derive the planar approximation expressions for the STO 
oscillation periods and compare them with the numeric 
results obtained without approximations. 



II. MODEL 

We will consider the case of a spin transfer (Fig. [T]) 
device with two magnetic layers, both described as 
macrospins. One of the layers has a fixed magnetic mo- 



ment with a direction given by a unit vector s. This layer 
acts as a spin polarizer. The other layers's magnetization 
is free to move and can be described by a macrospin mag- 
netic moment M(t) = Msn{t) where Mg is the constant 
saturation magnetization and n is a unit vector. 

The magnetic dynamics of the free layer is governed 
by the LLG equation with spin transfer term (see, e.g., 
Ref.fiF 
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where E is the total magnetic energy of the free layer, 
7 is the (positive) gyromagnetic ratio, a is the Gilbert 
damping, and the spin-transfer magnitude is given by 
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where V is the free layer volume, / is the electric current, 
and e is the electron charge. The factor I/e is positive 
when electrons flow into the free layer. Generally the 
spin polarization factor g{P, {n- s)) depends on the angle 
between the polarizer and the free layer and on the degree 
of spin polarization P.— In this paper we will employ 
the frequently used approximation g = const. In this case 
u is simply a rescaled current value. 

We consider the free layer with an easy plane 
anisotropy in the (x, y) plane and additional easy axis 
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FIG. 1: Spin transfer device with a fixed and free magnetic 
layers 



2 



anisotropy along the x direction. External magnetic field 
H is also applied along x. The magnetic energy is given 

by 
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where Kp and Ka are the easy axis and easy plane 
anisotropy constants. To simplify notation we will use 
the rescaled energy function 
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where the new constants ujp = ^Kp/M^, uja = ^Ka/Mg^ 
and h = jH have the dimensions of frequency. 

The fixed layer magnetization is assumed to be point- 
ing along the easy axis as well, s — +x. 



III. EFFECTIVE PLANAR DESCRIPTION 



As shown in Ref. |24| . when the inequalities ujp 3> uja 
and ujp ^ h hold, vector n moves close to the {x,y) 
plane and its magnetic dynamics can be described by an 
effective planar equation governing the behavior of its in- 
plane (azimuthal) angle 4>{t), measured from the a:-axis. 
The equation reads 
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(4) 



One can observe that it has a form of a Newton equation 
of motion for an "effective particle" with a mass l/ujp 
moving in the external one-dimensional potential Sef / (0) 
with a variable friction coefficient aeff{(j)). The effective 
energy and friction are found^** to be given by 
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The analogy between the effective planar equation and 
the Newton equation for a one-dimensional particle of- 
ten provides a qualitative understanding of the system's 
behavior. 

In the absence of spin transfer {u — 0) one has ae// — 
a > 0, and the planar equation predicts that the solu- 
tion (j){t) will approach a minimum of £eff{4>)- Indeed, 
a particle subjected to a viscous friction force will even- 
tually come to rest in one of the energy minima. This 
is guaranteed by a classical mechanics theorem, stating 
that the total energy 



etot{t) = ^ 



(6) 



is a decreasing function of time. The theorem is based 
on the relationship 



(7) 



and holds for Uef / > (case of conventional friction) . 

When the current is turned on, the effective friction 
can become negative for some values of cj). In the areas of 
negative friction the energy of the effective particle may 
increase, the theorem about the decrease of Etot {t) breaks 
down, and the particle motion does not necessarily have 
to stop in the local energy minimum. As a result, the 
effective particle can perform an oscillating motion which 
corresponds to the persistent precession state of the free 
layer. 



A. Description of oscillation regimes 

In the case of effective energy and friction give n by 
Eqs. [5] the effect of m 7^ was analyzed in Ref. 12J. We 
start with a review of these results. Denoting w = Wa + 
M^/wp, one finds that the effective energy has a minimum 
at = for /i > —Cj and cf) = tt ior h < uj. Both 
minima are present for \h\ < uj. For small currents the 
effective damping is still positive for all values of (j>. As 
the current increases, ae/ / first becomes negative at one 
of the minimum points = or = tt at the critical 
current values u = Tt^ij where 
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Negative a^f / at a point of energy minimum leads to the 
development of oscillations around this formerly stable 
equilibrium. As the current grows, the region of negative 
(Xeff becomes larger and the oscillations amplitude in- 
creases. The evolution of oscillations at large amplitudes 
proceeds differently for \h\ < uj and \h\ > uj. 

In the large field case, \h\ > Q, the energy has only 
one minimum. For dcfiniteness, consider the case of pos- 



itive h > UJ, with the minimum of e 



at 



and 



a maximum at (j) — tt (Fig. [2]). The effective damp- 
ing at the minimum point changes sign to negative at 
u ~ — ui. After that the system develops small oscilla- 
tions around (f> = (Fig.O^). As the current is made even 
more negative, the amplitude of these oscillations grows, 
until at a second critical current u = —U2{h) it becomes 
so large that the effective particle reaches the point of 
energy maximum (p = n. After that the particle starts to 
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FIG. 2: Oscillation states at /i > d). The solid line gives 
the effective energy eeff{4>)- Dashed lines with arrows show 
the motion of the particle in the potential profile. Horizontal 
dotted line gives the level of the total energy etot- (A) small 
oscillations at —U2 < u < —ui. (B) full rotations at u < —112. 



3 



(A) e 


(B) ^ 


?tot 




711 


1 1 V 


jj ^tot V 

// 1 1 , II Y 




,/ <!>,„ It -It -<]>,„ Vi 
8t 


,7 It (|) 



(C) 




FIG. 3: Oscillation states and switching at < /i < tl). The 
solid line gives the effective energy esff{<p)- Dashed lines with 
arrows show the motion of the particle in the potential profile. 
Horizontal dotted line gives the level of the total energy etot- 
(A) small oscillations at —U2 < u < —ui; (B) Switching to the 
(j) = n local minimum at —us < u < — U2; (C) Full rotations 
at u < — U3. 



make full circles in the easy plane (Fig. [2j3). The speed 
of the particle performing full circles becomes larger and 
larger as the current is decreased further below — U2. 

In the small field case, \h\ < uj, there are two minima 
of Sef / at = and </> = tt with a maximum point at the 
angle (j)m{h) between them (Fig. [3]). Considering again 
the case of positive field < h < u one finds the following 
picture. The (/) — minimum is destabilized at a nega- 
tive current u = —ui and small oscillations around this 
equilibrium are developed (Fig. . As the current is de- 
creased further, the amplitude of the oscillations grows 
until they reach the point of energy maximum at the sec- 
ond critical current u — —U2{h). For larger amplitudes 
the effective particle moves over the energy maximum 
into the basin of the (j) — t: minimum (Fig. [5}3). In this 
basin the effective damping is positive and the particle 
gradually looses its total energy, ending up in the = tt 
minimum. 

There is, however, yet another transition at the third 
critical current u ~ ~U3{h). Beyond the third thresh- 
old the particle reaches the (j) = 0m maximum with a 
velocity that is sufficiently high to allow it to move all 
the way to the next energy maximum at = 27r — (f>m, 
overcoming the friction force that attempts to stop it 
(Fig. In other words, the energy obtained in the 

negative damping area around = is sufficient to push 
the particle around the full circle from — 0„i to 27r — 0,„. 
At u < —U3{h) the particle can be either in the oscil- 
lation regime with full circle rotations, or at rest in the 
(j) = TT energy minimum. The actual state is determined 
by the history of the system. 

So far we have discussed small oscillations near the 
= equilibrium and their transformation into the full 
rotation regime. Positive current u > can destabilize 
the (j) = n equilibrium. Due to the symmetry of the 




FIG. 4: Switching diagram. The dashed and dash-dotted lines 
give the stability boundaries for the = tt, antiparallel (AP) 
and = 0, parallel (P) states respectively. These static states 
coexist with each other in the central rectangle of the diagram. 
The oscillatory regions are shown as grey areas, marked by 
the type of oscillation: small oscillations (SO) near P or AP 
points or full rotation regime (FR). Oscillations sometimes 
coexist with P or AP states. The diagram is symmetric in 
accord with the condition ([8]). 



problem, the critical currents Ui of the = tt equilibrium 
are related to those already found for (/) = as 



Ui{TT, h) — ~Ui{0, —h) i = 1, 2, 3 



(8) 



The switching diagram^* of our spin transfer system is 
shown in Fig. 01 It consists of the two 90-degree wedges 
of stability of the = equilibrium ("parallel state"), 
and (j) — TT equilibrium ("antiparallel state"). In addi- 
tion to the stability areas of the fixed points there are 
stability areas of the oscillating solutions (grey areas in 
Fig. U). The oscillating solutions of the planar equation 
are in a one-to-one correspondence with the free layer 
precession states previously found numerically by solv- 
ing the LLG equationi^i^ The small oscillations regime 
corresponds to the case of in-plane precession of vector 
n, and the full rotation regime corresponds to the out- 
of-plane precessioni^ 

As the current is increased deeper into the full rota- 
tion regime the deviations of n from the (x, y) easy plane 
plane grow, and the planar approximation becomes inap- 
plicable. The orbits of persistent precession evolv o^^'^^'^^ 
into small circles around the maximum of the magnetic 
energy ([3]) far away from the easy plane. Eventually this 
point is stabilized by spin transfer<^i^ 
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B. Small damping approximation 

Consider first the planar equation at u = 0. In 
the limit of a — >■ the friction term can be viewed as 
a perturbation on top of the frictionless motion. When 
friction is completely absent, the total energy is exactly 
conserved 



stotit) = etot(O) 



(9) 



For small friction Stot is conserved approximately, i.e., its 
relative change during one period of oscillations around 
the potential minimum is small. On a time scale of sev- 
eral periods one can approximately use equation (|9]). For 
longer time intervals one has to take into account that 
£tot{t) is a slowly changing function of time. 

In the case of u 7^ the same picture holds as long as 
the absolute value of ae// is small. The only difference 
is that now the total energy can either decrease or grow 
during one period of oscillations. Here we want to note 
that the smallness of |ae//| is already guaranteed by the 
smallness of a in the following sense. As discussed above, 
the oscillating regimes of spin transfer systems appear 
when the effective damping becomes negative. At the 
onset of the oscillation regime the Gilbert damping term 
and the spin transfer term in the expression for aef / are 
of the same order. That means that |ae/ / 1 and a have the 
same order of magnitude, unless the current is increased 
far above the threshold. 

The slow time change of etot can be calculated using 
the following procedure. We use Eq. ([7]) to find the en- 
ergy change during one period of oscillations 



A = - / aef/rdt . 
Jo 

Applying the approximate energy conservation during 
the period one can express the particle velocity as 

(j) « ±yj 2ujp{etot - ee//(0)) = /(</>, etot) (10) 

where etot is assumed to be constant, for example taken 
as the value of the total energy at the beginning of the 
period. The plus or minus sign in front of the expression 
is chosen according to the direction of particle motion. 
The electric current parameter u enters aef / as a linear 
correction, and Sgf / as a quadratic correction. It will be 
shown in Sec. IIII Dl that the correction to the effective 
energy can be dropped in our approximation and £e// 
can be substituted by s{(j)) = — (1/2) uj a cos^ (j) — hcoscj) 
in Eq. (|10p. i.e., one can use uja instead of a). 

Using approximation (|10p one finds both A and the 
oscillation period T as 



(^efs{<f))f{(t)) dcj) 



(11) 



T = dt 
Jo 



In these formulae the integrals over (/> are taken along the 
closed trajectory corresponding to one period of oscilla- 
tions in the absence of damping and spin transfer. The 
value of Etot determines both the integration limits, and 
the integrand. Consequently, equations ([TT|) give both A 
and T as the functions of etot- The slow time evolution of 
the total energy is described by an approximate equation 



dst 



dt 
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In the regime of persistent oscillations the total energy 
is constant, which implies A{etot) = 0. In Ref. [13 this 
condition was used to find the critical currents Ui by the 
following argument. The critical trajectories correspond- 
ing to each Ui were already described in Sec. IIII Al On 
those trajectories the value of etot is given by the value 
of Ee// at the turning points where = 0. Now the in- 
tegral ([TT]) can be calculated as a function of u, and the 
equation A(w) — gives a current threshold. 



C. Oscillation periods 



Here we use p2|) to calculate the period of oscillations 
with arbitrary amplitude. Our goal is to find the function 
T{u). The equation detot/dt = or A = holds for 
persistent oscillation regimes at any current, not just for 
the critical current values. At a given u this equation 
determines the value of etot and the endpoints of the 
integration contour in the first equation of the system 
(ITlT) . Knowing them, we can calculate T from the second 
equation in (|lip and find the function T[u). 

We will do the calculations for u < 0. According to the 
symmetry of the switching diagram (Fig. |4]) the u > 
results can be easily obtained from those for -u < 0. 



1. Case of h> uja 

Here one observes small oscillations for —U2 < u < —ui 
and full rotations for u < —U2- 

In the small oscillations regime the particle moves be- 
tween the points of maximum deviation ±0a (Fig. EK). 
The total energy is given by etot = £(0a) = 
(l/2)a;a cos^ 0a — hcos(j)a, and the expression (fTO|) spe- 
cializes to 
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Condition A = together with ([TT|) gives 
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R^cf), TT, h) as 



2u 
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with 



r4>a 



We can now express the current as a function of the os- 
cillation amplitude 4>a 



u{(j)a) 



2 K2Ua) 



(13) 



The second equation of (|lip gives the oscillation period, 
also expressing it as a function of (pa 



2K. 
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Using equations ([T^ and ([T^ one can make a parametric 
plot T(u) by varying (j)a from zero to tt. The integrals if^ 
can be expressed through special functions or, in practice, 
calculated by a computer algebra system. 

The boundary between the small oscillations and the 
full rotations is found at the current value 

- U2{h) ^ u{tt) ^ TT^-rr-, jT [h > Wa) , (15) 

2 A2(7r, ft) 

where we have explicitly indicated the dependence of the 
integrals on h. Note that this formula corrects formula 
(8) from Ref. Ill which is valid only for the values of h 
slightly above Wa, i.e., for h — Wa ^a- 

Below the —U2 threshold the oscillations happen with 
a full rotation of the angle (/). There are no endpoints 
here (Fig. [2j3) and we will parameterize the trajectory 
by the excess energy Ae = etot — £(7r) > 0. Equation 
(|T0| specializes to 



= f]^ i?(0, TT, h) + ^^^^ = V'i?(<^, Tr,h)+6 

where we have introduced S — 2iOpAe/fl. Repeating the 
steps that led to equations (IT3| and p4)) we get an anal- 
ogous pair of equations for the full rotation regime 



u{6) = 



aujp Li{S) 
2~L2{5) 



(16) 
(17) 



where the integrals Li are expressed through R 



Li{5) 



R{(f>) +5 d(j) 



L2{6) = / cos(/.yi?((/)) +(5 

J — TT 



Using equations (1161) and (1171) one can make a parametric 
plot T{u) by varying 5 from zero up. 



2. Case of \h\ < uja 

For small oscillations (—1*2 < u < — wi), the calcula- 
tions turn out to be identical with those performed in the 
previous section. Expressions ([T5)) and (ITi|) can be used 
to make a parametric plot T{u). The only difference is 
that now the oscillation amplitude 4>a changes form zero 
to the angle (pm of the energy maximum. 

The critical current ^2(^1) is found by setting (pa — 4>m- 
The position of the energy maximum is found from the 
equation ds/dp = which gives 



COS( 



h 



At the same time the expression for R can be rewritten 
as 



R = [ cos H 



COS H 



As a result, for pa — 4'm one is able to write an explicit 
formula for the square root 



\/R{<P, pm,h) = I cos -I- h/u!a\ 



(18) 



The integrals Ki^2 can be then taken and one gets the 
expression for the critical current^ 

sinPra + {h/uJa)p.m ^ . 

U2=aujp- — — - — , {\h\<uJa). (19) 

The third critical current u = —uslh) corresponds to 
the trajectory shown in Fig. [3p. The angle changes from 
—pm to 27r — p„ii and the total energy is equal to e(27r — 
p>m) because at the critical current the particle starts at 
—prn with zero velocity, makes a full rotation, and reaches 
27r — p„i with zero velocity. Since e{2Tr — p>m) = e{p>m), 
we find that the condition A = specializes to 



^"-'f'- ( 2ucos( 

0L^ 

-0™ V '^P 



\J R{(p, (j)m,h) dp ^ 



(the integrand is the same as in the case of small oscil- 
lations but the integration limits are different). Using 
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FIG. 5: Large field case. Comparison between the no- 
approximation numeric results (dots) and planar approxi- 
mation expressions (lines) for the oscillations periods. Here 
LOp = 100 LOa, h — 1.4 cja, ct = 0.01. The oscillation period 



diverges at u = — U2 

equation ([18]) we then find^i 
+ (/i/a;J(0 



U3 



n/2) 



tt/2) + (h/uja) sin0„ 



auj„ . , — . ', \ . ' , ' , (\h\<uja)- 

(20) 

For the currents exceeding the third threshold, u < 
—U3, there are no turning points and, similar to the 
case of large fields, we parameterize the trajectory by the 
scaled energy excess over the maximum: S = 2wp(e(0) — 
£(</'m))/^ > 0. The resulting parametric expressions 
for the current and period turn out to be identical to 
Eqs. and fTT]) , except that the meaning of R changes 
to i? = R{4>, (j>m, h) = {cos(j) + h/waY the definitions 
of the integrals i 1,2,3 • 

Overall, equation pairs (IT3l [T4)) and p6l [T7| give the 
planar approximation formulae for the periods of all spin 
transfer oscillations possible in our system. For small os- 
cillations regime with current just above the first thresh- 
old the results (fT4)) and p7)) give the frequencies converg- 
ing to y^(ZJ^r+"^iVpj to the Kittel's formula with a 
substitution ujp -\- h ^ ujp m. accord with our approxima- 
tion h <^ ujp. For other current values we will present 
the results of the theory as graphs (Figs. [S] and IH]) for the 
periods T(u, h) rather than for the frequencies since os- 
cillations periods are more directly interpreted in terms 
of effective particle analogy. 



D. Condition of small damping 

The validity of the small damping approximation for 
Eq. requires 



A < eto 



(21) 



We can make a rough estimate by noting that in the 
regime oi h ^ uja the value of R is of the order of one. 
The current magnitude is assumed to be such that the 
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FIG. 6: (Color online) Small field case. Comparison between 
the no-approximation numeric results (dots) and planar ap- 
proximation expressions (line) for the oscillations periods for 
ujp = 100 uja, h — 0.5 LJa- The oscillation period diverges at 
u — —U2 and it = —u^ with no oscillations between the two 
thresholds. Top panel (A) gives results for a = 0.001, where 
the agreement is very good. The bottom panel (B) shows a 
comparison for a — 0.01. Planar approximation works well 
for small precession, but the full rotations regime requires 
currents which are too large for the small friction approxima- 
tion (|16l I17p to be accurate. The exact solution of the planar 
equation (red (grey) dots) is still in very good agreement with 
the numeric solution of the LLG equation (black dots). 



spin transfer term in cce// is of the same order as bare a. 
Then the formula pT|) gives 

As for the total energy, we assume that it has the same 
order of magnitude as the magnetic energy 

£tot ^ £ ^ OJa ■ 

With the above estimates the small damping condition 
(pij) transforms into 



a < 



(22) 



This is a well known condition of small Gilbert damping 
in planar systems. 

We can now show that the current-dependent term in 
uj = u!a + u'^ /ujp is a small correction to uJa whenever the 
conditions ([2^ holds. Indeed, for currents of the same 
order as the critical current we have u ~ awp, so 



u 

Wo 



a UJp UJa 
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This estimate justifies our statement in Sec. IIIIBI 



IV. COMPARISON WITH THE 
NO-APPROXIMATION NUMERIC RESULTS 

Oscillation periods can be found numerically by solv- 
ing the LLG equation without approximations. In prac- 
tice this means rewriting the vector equation ^ in polar 
angles 9 and of the unit vector n = {nx,ny,nz} — 
{sin d cos 0, sin Q sin cos ff\ . One obtains a system 

; -f aFe 



(1 + a^) sin 6* 
Fe - aF^ 



with 

F^ 

Fn 



1 



- sin cos 0{ijjp -f LjJa cos <f>) — h cos 6 cos 
-usin0 , 

-Wa sin 9 sin cos (p — h sin 4> + u cos cos < 



(23) 



The system is solved numerically with Mathematica. 
When an oscillating solution exists, the motion of n ap- 
proaches the same precession cycle regardless of the ini- 
tial conditions, as long as they are its basin of attraction. 
The period of the precession motion is measured after 
the stationary regime is achieved. 

Fig. [5] compares the numeric LLG results with the pla- 
nar approximation formula for the high field case, h > LUa- 
One observes a very good correspondence. Near the crit- 
ical current u — —U2 the periods of oscillations become 
infinite and their frequencies drop to zerot^ This is easy 
to understand from the effective particle analogy: at the 
critical current the particle travels between the two max- 
ima of the energy profile that have the same height with 
the particle energy being equal to the potential energy at 
the maximum point. As it usually happens in such cases, 
the motion near the maximum is infinitely slow and the 
period is infinite. 

In the low field case, < h < oja, the comparison of 
numeric and approximate results is shown in Fig. [S] In 
the top panel (Fig. [6j'V) the Gilbert damping is set to 
a = 0.001, and the condition |Q;e//| <C \J L/JajijJp is well 
satisfied for all current magnitudes on the graph. The 
correspondence between the numeric LLG results and the 
small damping approximation to the planar equation is 
very good. As expected, the oscillations period diverges 
at It = —ui and u = — U3, with no oscillations between 
the two thresholds. 

Fig. |BJ3 shows results for a = 0.01. We observe good 
correspondence between the LLG numeric results and our 



theory for the small oscillations, but the full rotations 
regime shows appreciable differences which become very 
large near the critical current. Their origin is the break- 
down of the small damping approximation: for a — 0.01 
the currents in the full rotation regime are so large that 
the strong inequality |ae//| ^ \/i^a/^p is not well satis- 
fied. To prove that the small damping approximation is 
the source of the discrepancy we have solved the planar 
equation (jj]) numerically. The results (red (gray) dots in 
Fig. |B|3) correspond very well to those obtained directly 
from LLG (black dots in Fig. [Bp). 

It is also instructive to compare the cases of large 
and small fields with the same value of Gilbert damp- 
ing a — 0.01 (Fig. [S] and Fig. HP). One can see that the 
validity region of the small damping approximation ex- 
tends at least to m ~ — 2 in the case of large field. At the 
same time in the case of small fields this approximation is 
visibly violated for m « — 2. The relative fragility of the 
small damping approximation in the full rotation regime 
in small fields can be traced to the presence of two en- 
ergy maxima, (j) = ±(j)mi instead of just one maximum, 
= TT, in the large fields. The situation calls for more 
work on the approximate solutions of the effective planar 
equation Q with variable damping. 



V. CONCLUSIONS 

We have shown that the planar approximatio n^'^i^^ 
gives good results for the frequencies of spin transfer os- 
cillators with dominating easy plane anisotropy. Analytic 
expressions for the oscillation periods were derived in the 
limit of small Gilbert damping. The mechanical analogy, 
associated with the effective planar equation, provides a 
qualitative understanding of different precession regimes 
and naturally explains the singular behavior of the pre- 
cession frequency near the transition between in-plane 
and out-of-plane precessions. 

An important advantage of the planar approximation 
is the fact that the problem of finding the unperturbed 
trajectory is drastically simplified. All one-dimensional 
trajectories are straight lines completely characterized by 
their endpoints. This property should help to develop 
theories of the large- angle precession regimes of planar 
spin torque oscillators in the presence of temperature 
fluctuations or other noise sources. 
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